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Abstract 

Different interoceptive systems must be integrated to ensure tliat multiple homeostatic insults evol<e appropriate 
behavioral and physiological responses. Little is known about how this is achieved. Using C. elegans, we dissect cross- 
modulation between systems that monitor temperature, O2 and CO2. CO2 is less aversive to animals acclimated to 1 5°C than 
those grown at 22°C. This difference requires the AFD neurons, which respond to both temperature and CO2 changes. CO2 
evokes distinct AFD Ca^^ responses in animals acclimated at 15 'C or 22°C. Mutants defective in synaptic transmission can 
reprogram AFD CO2 responses according to temperature experience, suggesting reprogramming occurs cell autonomously. 
AFD is exquisitely sensitive to CO2. Surprisingly, gradients of 0.01% C02/second evoke very different Ca^^ responses from 
gradients of 0.04% C02/second. Ambient O2 provides further contextual modulation of CO2 avoidance. At 21% O2 tonic 
signalling from the 02-sensing neuron URX inhibits CO2 avoidance. This inhibition can be graded according to O2 levels. In a 
natural wild isolate, a switch from 21% to 19% O2 is sufficient to convert CO2 from a neutral to an aversive cue. This sharp 
tuning is conferred partly by the neuroglobin GLB-5. The modulatory effects of O2 on CO2 avoidance involve the RIA 
interneurons, which are post-synaptic to URX and exhibit C02-evoked Ca^"^ responses. Ambient O2 and acclimation 
temperature act combinatorially to modulate CO2 responsiveness. Our work highlights the integrated architecture of 
homeostatic responses in C. elegans. 
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Introduction 

To maintain a constant internal milieu animals use internal 
sensory receptors to monitor cues sucti as COj/pH [1], [2], 
temperature [3], and osmolality [4], These interoceptors counter 
changes in internal mUieu by coordinating homeostatic responses 
that alter physiology and behavior [5]. Cross-talk between 
different interoceptive systems is likely to be important to ensure 
an integrated homeostatic response by the animal to multiple 
homeostatic insults. However, relatively little is known, at the 
molecular and circuitry levels, about how such cross-talk is 
encoded. 

In vertebrates electrophysiological studies have identified cell 
populations and circuits that respond to homeostatic imbalance in 
O2, CO2/PH and temperature. The neurons comprising these 
circuits are only beginning to be resolved, and the molecular 
mechanisms controlling their responses are poorly understood. 
Nevertheless, studies in several animals suggest that cross- 
modulation of homeostatic responses is important for survival. In 
panting mammals, a rise in core body temperature elicits increased 
ventilation rate to help cooling, even though this causes temporary 
alkalosis of the blood due to excessive blowing off of CO2. This 
over-ride appears to be achieved by changing the set-point at 
which CO2 sensors inhibit ventilation when [CO2] decreases, but 



the mechanisms involved are unclear [6]. In the mouse, recent 
work has shown that suppressing the activity of serotonergic 
neurons impairs both respiratory and body temperature control, 
although whether the same or different sub-populations of neurons 
mediate these effects is unclear [7,8]. In mammals, the drive to 
increase ventilation rate is stimulated more strongly when animals 
simultaneously experience a drop in O2 and a rise in CO2 [9]. 

In invertebrates, such as the free-living nematode C. elegans, 
behavioral mechanisms that counter homeostatic imbalance are 
particularly important, since the animal's buffering capacity is 
limited. C. elegans responds to variation in temperature, O2 and 
CO2 by mounting sophisticated behavioral responses. Exposure to 
temperatures above or below the range in which C. elegans can 
grow elicits strong avoidance responses [10]. When navigating 
thermal clines in which it can thrive, ~15°C to 25°C, C. elegans 
migrates to the temperature at which it grew recently, as long as 
this was not associated with starvation [11,12]. These responses 
require the animal to memorize its recent temperature experience 
and to change this memory when temperature or nutrient 
conditions change. A neural circuit that subserves these behaviors 
has been identified, and involves the thermosensory neurons AFD 
and AWC [13-16]. Temperature experience alters the thermo- 
sensing properties of AFD neurons: in animals acclimated to 
higher temperatures, the threshold at which a temperature rise 
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Author Summary 

Many animals are either attracted or repelled by carbon 
dioxide. We show that the way C. elegans responds to CO2 
depends on the temperature it has acclimated to and the 
oxygen tensions it is experiencing. The effects of acclima- 
tion temperature are mediated by a temperature-sensing 
neuron called AFD that also responds to CO2. The 
responses evoked in AFD by a change in CO2 concentra- 
tion are reprogrammed by acclimation temperature. This 
reprogramming does not appear to require synaptic input 
from other neurons. O2 modulates CO2 avoidance by 
setting the activity of the tonically signalling O2 sensor 
URX. A switch from 21% to 19% O2 is sufficient to convert 
CO2 from a neutral stimulus to an aversive one in a C. 
elegans wild strain. Modulation of CO2 responses by O2 
cues requires the interneuron RIA which itself responds to 
changes in CO2 and is directly post-synaptic to URX. CO2 
gradients are likely to be common in rotting fruit where 
Caenorhabditis live. Such gradients could be associated 
with food, pathogens, conspecifics or predators of C. 
elegans. The value of CO2 as a sensory cue thus depends 
crucially on context. This may explain the remarkable 
complexity of CO2 sensing in C elegans. 



evokes a Ca response in AFD occurs at correspondingly higher 
temperatures [17,18]. This plasticity allows animals to respond 
homeostatically to external temperature fluctuations, by seeking 
and remaining at temperatures they are acclimated to. 

C. elegans also displays responses to variation in [O2] , and avoids 
both high and low O2 [19]. Wild-caught C. elegans strongly avoids 
2 1 % O2, both on and ofT food, and burrow to escape from the 
surface [20]. This avoidance response is sculpted by 02-sensing 
neurons in the body cavity called AQR, PQR and URX 
[20,21,22]. When [O2] levels rise towards 21% the AQR, PQR 
and URX neurons become activated, by a mechanism involving 
the atypical soluble guanylate cyclases GCY-35/GCY-36. The 
tuning of the O2 response is sharpened by a neuroglobin expressed 
in AQR, PQR and URX neurons, called GLB-5, that suppresses 
neuronal activity when ambient [O2] falls just below 21% [20,23]. 
The AQR, PQR and URX neurons are aU tonic receptors: they 
show sustained signalling as long as [O2] is high [24]. This tonic 
activity stimulates sustained rapid movement until animals 
encounter a preferred lower [O2] environment. 

C. elegans also avoids elevated CO2 [25,26]. As in vertebrates, 
high [CO2] is harmful to C. elegans, reducing brood size and 
disrupting muscle structure [27]. An array of sensory neurons 
mediates CO2 avoidance behavior [28]. This network includes the 
temperature sensor AFD, the major gustatory neuron ASE, and 
the BAG neurons, which are also activated by decreasing O9 levels 
[22]. 

Here we investigate how the temperature and O2 sensing 
systems of C. elegans modulate the distributed circuit that mediates 
responses to CO2. 

Results 

Previous temperature experience sets CO2 avoidance in 
C. elegans 

To examine if temperature can modify C. elegans' responses to 
CO2 we grew N2(Bristol) animals at 22°C and compared their 
behavior in CO2 gradients at 15°C and 22°C (Figure lA, B) 
[25,28]. CO2 avoidance at the two temperatures was similar when 
animals navigated 3%-0% and 5%-0% CO2 gradients. However, 
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Figure 1. CO2 avoidance is modulated by acclimation temper- 
ature. A. Assay for C. elegans CO2 responses. Animals navigate a 
defined CO2 gradient in a microfluidic device. The chemotaxis index is 
calculated by counting animals in two halves of the device, using the 
formula shown. B-D. Chemotaxis indices for animals cultivated at either 
15°C or 22 C and assayed in different CO2 gradients at either 15"C or 
22°C. **, p<0.01; n.s., not significant, Student's f-test. E. A mutation in 
ffx-!, which is specifically required to confer AFD neural identity, 
disrupts modulation of CO2 avoidance by acclimation temperature. 
Assays were performed in 3%-0% CO2 gradients. **, p<0.01; n.s., not 
significant, Student's f-test. 
doi:1 0.1 371/journal.pgen.1 004011 .gOOl 



animals in a 1-0% CO2 gradient avoided the high CO2 half of the 
microfluidic device more strongly when assayed at 15°C compared 
to 22°C (Figure lA). 

C. elegans can retune its temperature preference according to the 
temperature to which it is acclimated [13,29]. This behavior is 
encoded in AFD [17,18], a neuron that also responds to CO2 [28]. 
We therefore examined how previous temperature experience 
altered subsequent CO2 responses. We grew animals at 15°C or 
22°C, and assayed their CO2 responses at each temperature. 
Strikingly, previous temperature experience altered CO2 avoid- 
ance. Animals grown at 15°C avoided CO2 less strongly than 
animals grown at 22°C, regardless of whether the assay 
temperature was 15°C or 22°C (Figure IB— D). Animals grown 
at 15°C showed weaker CO2 avoidance even when exposed to 
relatively high CO2 levels, 5% (Figure IB-D). Thus, the 
temperature to which C. elegans has acclimated helps determine 
the aversiveness of CO2. 
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Acclimation temperature does not reprogram CO2 

responses in AFD-defective mutants 

We investigated if the AFD neurons helped to reprogram CO2 
avoidance behavior according to acclimation temperature. The 
ttx-1 (diermo&«s defective) gene encodes a member of the OTD/ 
OTX subclass of homeodomain transcription factors [30]. 
Mutation.s in ttx-1 selectively disrupt AFD specification, and 
confer a thermotaxis-defective phenotype. Loss of ttx-1 also 
reduces CO2 avoidance in animals navigating COg spatial 
gradients [28]. If AFD neurons were important for temperature 
regulation of COj avoidance responses, then ttx-1 mutants would 
display similar CO2 avoidance regardless of cultivation tempera- 
ture. As shown previously, ttx-1 mutants grown at 22°C only 
avoided CO2 weakly [28], resembling wild-type animals grown at 
15°C (Figure IE). This defect was rescued by a wild-type ttx-1 
transgene (Figure IE). By contrast, loss of ttx-1 did not alter the 
COj-avoidance behavior of animals cultivated at 15°C. These 
data suggest AFD is required for acclimation temperature to 
modify COj aversive responses. 

Acclimation temperature re-programs the CO2 

responsiveness of AFD 

Acclimation temperature sets the response threshold of AFD 
neurons to warming [17]. This prompted us to investigate whether 
acclimation temperature also alters the CO2 responsiveness of 
AFD. To measure C02-evoked Ca^"*" responses in AFD we 
expressed the genetically encoded Ca^"*" sensor cameleon YC3.60 
[31] from the gcy-8 promoter [32]. For our recordings we used 
animals acclimated to 15' C or 22 C, but maintained animals at 
22°C while we imaged them. In animals acclimated to 22°C high 
CO2 evoked in AFD the complex Ca^^ response described 
previously (Figure 2A) [28] . This typically consisted of an initial 
slight drop in Ca^""" when CO2 levels rose, followed by a rise in 
Ca^^ to above pre-stimulus levels, and finally, when the CO2 
stimulus was removed, a Ca^^ spike that rapidly decayed back to 
baseline. By contrast, animals acclimated to 15°C exhibited a 
simple response: a rise in Ca^"^ when CO2 levels rose, and a fall 
when CO2 was removed (Figure 2B). These data suggest that the 
previous temperature experience of C. ekgans reconfigures the CO2 
response properties of AFD neurons. 

To investigate if this retuning was driven by the intrinsic 
temperature-sensing properties of AFD neurons, or required pre- 
synaptic input, we imaged the Ca^^ responses of AFD neurons to 
CO2 in snh-1 (synaptoArevin-i) mutants, which are defective in 
synaptic transmission [33]. C02-evoked responses in AFD 
neurons were not altered in snb-1 animals compared to wild type, 
regardless of acclimation temperature (Figure 2C, D). These data 
suggest that the temperature experience can retune the CO2 
response properties of AFD neurons when synaptic signalling is 
defective. 

We characterized the response properties of the AFD neurons 
further. Previously, we had only exposed animals to sharp changes 
in CO2 that occurred within 1-2 s, and we always returned 
animals to 0% CO2 between stimuli [29]. To examine AFD 
responses to rises in CO2 from non-zero levels, we subjected 
animals acclimated to 22°C to a stimulus train involving multiple 
CO2 switches, namely 0%-l %-3%-5%-3%-l %-0%. Whenever 
CO2 levels increased, we observed an initial drop in Ca^""" followed 
by a rise in Ca^"^ (Figure 2E). Whenever CO2 levels decreased, we 
observed a spike of Ca^"*" that rapidly returned to baseline. This 
pattern of CO2 evoked Ca^""" response suggests that AFD can 
encode whether an animal is moving towards higher or lower 
CO2. 



Previous work has identified one potential molecular sensor for 
CO2, the transmembrane guanylate cyclase gcji-9 [34]. We 
compared C02-evoked responses in AFD neurons in wild type 
and giy-9 mutants. We observed no dilference in the response, 
suggesting that molecules other than GCY-9 confer CO2- 
responsiveness to AFD neurons (Figure SI). 

AFD responses to CO2 are reconfigured by the steepness 
of the CO2 gradient 

The ubiquity of CO2 suggests that its value as a cue is likely to 
depend not only on context (such as temperature) but also on the 
shape of the CO2 stimulus. Very rapid change in COj levels may 
convey a different meaning from a very gradual change. In our 
behavioral expcrimc'nts, animals navigated shallow CO2 gradients 
and encountered changes in the order of 0.01% CO2 per second 
(depending on speed and direction of travel in the gradient). To 
examine if AFD could respond to such shallow CO2 gradients, we 
exposed animals cultivated at 22°C to gradual linear increases and 
decreases in CO2 concentration at rates of 0.04% and 0.01% per 
second (Figure 3A,B). AFD responded to both these CO2 
gradients, but with very different response patterns. Gradients of 
0.04% C02/second evoked AFD Ca^"*" responses reminiscent of 
those elicited by sharp changes in CO2 (>1% C02/second; see 
Figure 2): Ca^^ levels decreased while CO2 was slowly rising to 
5%, then rose sharply as CO2 levels stabilized at 5%. When we 
gradually reduced CO2 levels back to 0%, Ca^^ levels spiked, 
returning to baseUne when animals were in 0% CO2 (Figure 3A). 
By contrast, gradients of 0.01% C02/second evoked a series of 
Ca^"^ spikes while CO2 levels were rising (Figure 3B). Ca^"*" levels 
tended to return to baseline when CO2 levels stopped rising, but 
spiking resumed when CO2 levels started falling. This spiking 
pattern disappeared when we imaged Ca^"^ responses evoked by 
the same 0.01% C02/second gradient in animals acclimated to 
15°C (Figure 3C). In these animals responses were more similar to 
those evoked by steeper CO2 gradients in animals acclimated to 
15°C (compare Figure 3C to Figure 2B). These data indicate that 
AFD neurons respond to both rapid and slow changes in CO2, but 
with different response patterns. The also highlight complexity in 
how AFD encodes CO2 stimuli. 

Ambient O2 levels regulate C. elegans CO2 avoidance 
behaviour 

To investigate further how different homeostatic responses are 
integrated, we examined if CO2 avoidance behavior was 
modulated by different background ambient [O2]. In body fluids 
and many ecological niches low [O2] coincides with high [CO2], 
and, conversely, 21% O2 is associated with low CO2. Cross-talk 
between the two gas sensing circuits could enable C. ekgans to 
recognize and respond appropriately to such environments. 

To examine this possibility, we placed N2 animals in 
microfluidic chambers containing gradients of CO2 at different 
fixed conc:entrations of O2. As expected, increasing [CO2] elicited 
increasing avoidance behavior: C. elegans avoidcxi j"'o CO2 more 
strongly than 3% or 1% CO2 (Figure 4A) [25,26]. Moreover, CO2 
avoidance was influenced by the background ambient O2 
concentration. N2 animals navigated down CO2 gradients more 
strongly when ambient O2 concentration was 11%, than when it 
was 21%. Increased avoidance was particularly striking when 
animals navigated shallow gradients of 1-0% CO2 (Figure 4A). 
Such shallow CO2 gradients are likely to be ecologically relevant 
in the rotting habitats where C. elegans thrives. 

To test the dynamic range of O2 regulation, we asked if 
increasing [O2] to above 21% could further suppress CO2 
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Figure 2. Acclimation temperature alters C02-evoked Ca responses in AFD neurons. In animals cultivated at 22 C a rise and fall in CO2 
evokes a complex Ca^* response in AFD neurons (A). Ca'^* initially falls when CO2 begins to increase, then rises. When CO2 levels fall, there is a Ca^* 
spike. By contrast, animals cultivated at 15 C show a simple response to the same stimulus (B). C-D The effect of acclimation temperature on CO2- 
evoked Ca'^* responses in AFD neurons is unaltered in snb-1 mutants defective in synaptobrevin. E. Ca^* responses evoked in AFD by a 0%-1%-3%- 
5%-3%-1%-0% CO2 stimulus train in animals acclimated to 22'C. Shading highlights switch times. Acclimation temperature is shown for each panel 
under genotype. 

doi:1 0.1 371/journal.pgen.1 00401 1 .g002 



avoidance. Although this is unphysiological, previous studies have 
shown that C. elegans can grow and reproduce in even 100% O2 
without any apparent adverse effects [35]. Since C. elegans only 
weakly avoided 1 % CO2 in 2 1 % Oj, we used a steeper 3-0% CO2 
gradient, to improve our dynamic range. Increasing ambient [O2] 
to 50% significantly suppressed avoidance of 3% CO2 (Figure 4B). 
These data suggest that ambient O2 concentration provides a 
contextual cue to modulate C. ekgans avoidance of CO2. 

Tonically signalling O2 sensors inhibit CO2 avoidance at 
high ambient [O2] 

Our results suggested that 02-sensing neurons or neuroendo- 
crine cells persistentiy signal O2 concentration to modify the 
activity of CO2 transducing circuits. Previous studies have shown 



that the AQR, PQR and URX O2 sensors signal tonically when 
ambient [O2] is close to 21%, and become progressively less active 
as [O2] falls [24]. The 02-evoked Ca^"^ responses of these neurons 
requires the atypical soluble guanylyl cyclases GGY-35 and GCY- 
36, which appear to be O2 sensors [20,22,36]. In gcy-35 or gcy-36 
loss-of-function mutants the Ca^""" levels in the O2 sensing neurons 
reported by cameleon YC3.60, are low, resembling those found in 
wild type animals kept at low [O2] [24] . To test if tonic signaUing 
by AQR, PQR and URX neurons persistently repressed CO2 
avoidance in high [O2], we compared the CO2 avoidance of wild 
type, gcy'35, and gcy-36 mutants at 2 1 % and 1 1 % O2. In 1 1 % O2 
gcy-35 and gcy-36 mutants avoided CO2 hke N2 controls. 
However, whereas increasing background O2 levels to 21% 
inhibited the CO2 avoidance behavior of wild type animals, it had 
no effect on gcy-35 or gcy-36 mutant animals (Figure 5A). These 
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Figure 3. Shallow and steep CO2 gradients evoke qualitatively different Ca'^* responses in AFD. A. Ca^^ responses evoked in AFD by CO2 
switches indicated at top, involving linear 0-5% and 5%-0% CO2 gradients occurring over 2 minutes. This corresponds to a rate of change of 0.04% 
C02/second. The upper part of the panel shows traces obtained from 10 randomly selected individual AFD neurons; an average trace is plotted at the 
bottom. Animals imaged in this panel were acclimated to 22 C. B, C. Ca^* responses evoked in AFD by CO2 switches indicated at top, involving linear 
switches from 0-5% and 5%-0% CO2 occurring over 8 minutes. This corresponds to a change of 0.01% C02/second. The upper part of the panels 
shows traces obtained from 10 randomly selected individual AFD neurons; average traces are plotted at the bottom. Animals imaged in (B) were 
acclimated to 22°C; those in (C) were acclimated at 15°C. For each panel, individual and average traces are at the same scale. The scale bar in each 
panel represents 0.4 YFP/CFP ratio unit. 
doi:1 0.1 371 /journal.pgen.1 00401 1 .g003 



data suggest that tonic signalling from one or more of the AQR, 
PQR and URX O2 sensors represses CO2 avoidance at high O2 
concentrations. 



To confirm our results, we rescued the gcy-3 6 mutant phenotype 
using cell-specific promoters. Expressing gcj-36 cDNA from its 
own upstream sequence, which drives expression in AQR, PQR 
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and URX, restored to gcy-36 mutants reduced CO2 avoidance at 
2 1 % O2 (Figure 5B), gcy-36 mutants expressing gcj-36 cDNA from 
the gcy-32 promoter, which also drives expression in AQR, PQR 
and URX, gave similar rescue (Figure 5B). Expressing gcy-36 
cDNA from the Jlp-8 promoter, which drives expression in URX 
(and AUA and PVM) neurons but not in AQR and PQR also 
rescued the 02-regulated CO2 avoidance phenotype of gcy-36 
mutants. These results suggest that tonic signalling by the URX 
02-sensing neuron can persistently suppress CO2 avoidance while 
O2 levels are high. 

To extend our results we also examined the consequence of 
deleting gcy-32 and gcy-34, atypical soluble guanylate cyclases 
expressed in AQR, PQR and URX neurons whose activities are 
also likely to be modulated by O2, but whose deletion only subtly 
alters 02-evoked behaviors. We observed no effects of these 
deletions on O2 regulation of CO2 avoidance (Figure S2). We did 
however observe a slight decrease in CO2 avoidance at 1 1 % O2 in 
mutants defective in gcy-33, an atypical soluble guanylate cyclase 
required for the BAG sensory neurons to respond to decreases in 
O2 levels (Figure S2) [22]. BAG is also a major CO2 sensor [28] 
[34]. 

The npr-1 and glb-5 genes modulate CO2 avoidance by 
O2 

O2 responses in the standard laboratory N2 strain differ from 
those of aggregating wild C. elegans, due to genetic differences that 
have evolved during domestication [19,20,23,36,37]. N2 animals 
harbor a gain-of-function allele of the npr-1 neuropeptide receptor 
that inhibits signalling output from 02-sensing circuits in feeding 
animals. N2 animals also carry a loss-of-function mutation in the 
neuroglobin that increases the excitability of the AQR, PQR 
and URX O2 sensors. 

We investigated if variation at npr-1 and glh-5 altered O2 
modulation of CO2 avoidance. In N2 animals, stepwise increases 
in O2 from 11% to 21% caused stepwise decreases in CO2 
avoidance (Figure 6A). Animals defective in both the npr-1 
receptor and the glb~5 neuroglobin (i.e. npr-1 mutants) were 
attracted to CO2 at 21% O2, but became progressively more 
repelled by CO2 as O2 concentrations fell. A functional glb- 
5 (Hawaii) allele made CO2 more aversive to npr-1 defective 
animals: decreasing [O2] still stimulated CO2 avoidance, but at 
each concentration tested glb-5; npr-1 animals avoided CO2 more 
strongly than npr-1 animals (Figure 6A). Adding the functional glb- 
5(Hawaii) allele to N2 animals bearing the npr-1 gain-of-function 
allele did not significandy change their CO2 avoidance behaviour 
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Figure 5. Disrupting gcy-35 or gcy-36 confers strong CO2 
avoidance regardless of ambient O2. A. gcy-35 or gcy-36 mutants 
strongly avoid the high CO2 half of a 1-0% CO2 gradient regardless of 
ambient O2. Statistics refer to comparisons to N2 at 21% O2. ***, 
p<0.001, Anova, Bonferroni corrected p-value. None of the strains apart 
from N2 show significant differences between assays carried out at 21% 
and 1 1 % O2 (Student's f-test). B. The C02-avoidance phenotype of gcy- 
36 mutants can be rescued by expressing gcy-36 cDNA in AQR, PQR and 
URX, using gcy-32 or gcy-36 promoters, or in URX alone, using the flp-8 
promoter. **, p<0.01, ***, p<0.001, Anova, Bonferroni corrected p- 
value. 

doi:1 0.1 371/journal.pgen.1 00401 l.gOOS 

at any O2 tensions. Thus, variation at the glb-5 and npr-1 genes, 
which alter O2 sensing circuits, changes the extent to which O2 
levels modifies CO2 aversiveness. 

To investigate how O2 modified CO2 avoidance in a non- 
domesticated C. elegans strain, we examined the responses of 
animals from the Hawaiian CB4856 isolate. As reported 
previously [23,25,26], the Hawaiian strain showed weaker CO2 
avoidance than N2 at 21% O2. Reducing O2 levels to 19% was 
sufficient to strongly stimulate CO2 avoidance in Hawaiian 
animals, and further decreases in [O2] had no significant effects 
(Figure 6A, C). Together, these data suggest that the Hawaiian 
animals do not avoid CO2 when O2 is at 21%, i.e. when animals 
are at the surface, and but that very small decreases in O2 are 
sufficient to increase C02-avoidance behavior. The sharp tuning 
of CB4856 responses to CO2 by O2 levels appears to involve the 
natural alleles of npr-1, npr-1 215F, the glb-5(Haw) alleles. 

To shed further light on the genetic control of this cross-talk of 
CO2 and O2 responses, we examined how knocking out the 
soluble guanylate cyclases ^g^-^J and ^tj-^ff altered CO2 responses 
in different genetic backgrounds. Knocking out either soluble 
guanylate cyclase strongly stimulated CO2 avoidance in npr-1 
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Figure 6. Re-configuring O2 sensing circuits by altering the npr- 
7 and glb-5 genes alters CO2 avoidance behavior. A. Tuning of 
CO2 avoidance behavior by different O2 concentrations in N2 (Bristol), 
CB4856 (Hawaiian), npr-Uad609), glb-5(Haw); npr-1 (ad609), and glb- 
5(Haw) animals. All assays used a 1-0% CO2 gradient. Statistical 
comparisons are to the N2 response at the corresponding O2 
concentration, *** p<0.001; ** p<0.01; *p<0.05 (Anova, p value 
protected Fisher's LSD). B. gcy-35 and gcy-36 mutants strongly avoid 
CO2 regardless of genotype at the npr-1 locus. Statistical comparisons 
are to the npr-1 response at the corresponding O2 concentration. ***, 
p<0.001, *♦, p<0.01, Anova, Bonferroni corrected p value). C. CB4856 
(Hawaii) animals defective in gcy-36 strongly avoid CO2 regardless of O2 
levels. Statistical comparisons are to the CB4856 response at the 
corresponding O2 concentration. ***, p<0.001, **, p<0.01, Anova, 
Bonferroni corrected p value). D, E. Tonic Ca^* levels in URX neurons of 
glb-5(Haw); npr-1 animals kept at 21% O2 and 17% O2 is lower than Ca'^* 
levels in URX in npr-1 animals kept at the corresponding O2 
concentrations. Ca^* measurements were made using cameleon YC2.60. 
doi:1 0.1 371 /journal.pgen.1 004011 .g006 

animals: the avoidance behaviour of gcj-35; npr-1 ot gcj-36; npr-1 
animals resembled that of gcy-35 or gcy-36 mutants, and of N2 
animals at 1 1 % O2 (Figure 6B), We also examined the effect of 
disrupting ^(y-j?^ in the Hawaiian genetic background (Figure 6C). 
CB4856 animals defective in gcy-36 avoided CO2 much more 
strongly than CB4856 controls, and changing ambient O2 had 
little effect on their CO2 responses (Figure 6C). Thus, the 
modulation we describe in domesticated N2 also occurs in wild 
aggregating C. elegans. Expressing cDNA encoding the npr-1 215V 
allele found in N2 animals in the AQR, PQR and URX neurons, 
using the gcy-32 promoter, restored N2-like behaviour to npr-1 



mutants (Figure 6B). Thus, npr-1 acts in the 02-sensing neurons 
themselves to counter the inhibitory effect of high O2 on CO2 
avoidance. 

To provide a neural explanation for why npr-1 animals avoided 
C02 less than glb-5(Haw); npr-1 animals at 17%, 19% and 21% O2 
(Figure 6A, /)<0.0001, Anova, Bonferroni-corrected p value at all 
three O2 values), we compared tonic Ca^""" signalling in URX at 
different O2 concentrations. While URX Ca^"*" levels were similar 
in npr-1 and glb-5; npr-1 animals at 7% O2, Ca^**" was higher in npr- 

1 than \n glb-5; npr~l animals at 21% and 17% O2, consistent with 
greater inhibition of CO2 avoidance by URX signalling at these 

02 concentrations (Figure 6D, E). 

O2 can modulate CO2 avoidance in animals defective in 
AFD and BAG CO2 sensors 

CO2 avoidance in C. elegans is mediated by a distributed set of 
sensory neurons that includes the BAG O2 sensor, the AFD 
temperature sensor, and the ASE gustatory neuron [28,34]. To 
examine if O2 levels modified C02-evoked Ca^^ responses in any 
of these neurons we imaged their responses at 1 1 % and 2 1 % O2 
concentrations using the YC3.60 sensor (Figure S3A-C). We did 
not observe any differences between C02-evoked responses at the 
two O2 concentrations in any of the three neurons under our 
imaging conditions. This suggests either that O2 modulation 
occurs downstream of these sensory neurons, or that our imaging 
conditions limit our ability to observe modulation by O2. 

O2 input could selectively modulate the CO2 responses 
mediated by one C02-sensing neuron, or it could modulate 
circuits involving multiple CO2 sensors. To examine these 
possibilities, we specifically disrupted AFD and/or BAG 
function in N2 animals, and measured CO2 avoidance at 21% 
and 11% O2. Genetically abating BAG neurons or disrupting 
AFD specification by mutating the ttx-1 transcription factor, or 
doing both, reduced CO2 avoidance at 11% O2, but did not 
abolish modulation by ambient O2 levels (Figure 7). These data 
suggest that O2 levels either modulate the output from several 
CO2 sensors, or exert their effects on unidentified CO2 sensors, 
or both. 

RIA intemeurons are part of the circuit mediating O2- 
modulated CO2 avoidance 

To dissect further how 02-sensing neurons modulated CO2 
responses, we sought mutations that disrupted O2 modulation 
without abrogating CO2 responsiveness. One such mutation we 
identified was ttx-7, which disrupts a m))0-inositol- 1 -monopho- 
sphatase [38]. ttx-7 mutants showed only mild defects in CO2 
avoidance when assayed at 21% O2 (Figure 8A-C). The 
chemotaxis index of ttx-7 mutants was not significantly different 
from that of N2 controls when animals were assayed in 1-0% and 
5-0% CO2 gradients; we only observed a small but significant 
decrease in CO2 avoidance when ttx-7 mutants were assayed in 3- 
0% CO2 gradients. However, ttx-7 mutant animals did not 
increase their CO2 avoidance when assayed at 1 1 % O2, regardless 
of the CO2 gradient we used (Figure 8A-C). ttx-7 mutants behaved 
indistinguishably from N2 animals when assayed in O2 gradients 
(Figure S4), suggesting they were not generally defective in O2- 
evoked responses. 

To confirm that the defect in 02-dependent modulation of CO2 
avoidance was due to the ttx-7 mutation, we showed we could 
restore strong CO2 avoidance at 11% O2 to ttx-7 mutants by 
expressing ttx-7 cDNA from the ttx-7 promoter (Figure 8D), 
Together, these data suggest that ttx-7 mutants can sense and 
respond to O2 but cannot communicate information about 
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Figure 7. Ambient O2 can modulate CO2 avoidance in animals 
lacking BAG and AFD CO2 sensors. Animals in which BAG neurons 
are ablated by specific expression of egl-1 caspase, and AFD neurons are 
defective due to loss of ttx-1, retain 02-modulation of CO2 avoidance. 
egl-1 expression in BAG neurons is driven by the flp-17 promoter. ^-^^ 
p<0.001, Student's f test, comparing a strain's responses at 21% and 
11% O2. *** p<0.001; ** p<0.01; * p<0.05, Anova, Bonferroni corrected 
p value, comparing responses to that of N2 at the same O2 
concentration. 

doi:1 0.1 371 /journal.pgen.1 00401 1.g007 



ambient [O2] to the appropriate circuits tiiat mediate CO2 
responses. 

To identify neurons where tlx- 7 acts to promote CO2 avoidance 
at low [Oj] we rescued the tlx- 7 COj avoidance phenotype by 
driving ttx-7 cDNA in small subsets of neurons. We focussed on 
neurons that receive synaptic input from the URX O2 sensors, 
since our gcy-36 rescue experiments imphed that URX was 
sufficient for O2 to modulate CO2 avoidance (Figure 5B). URX 
neurons make several synapses onto the RIA interneurons [39]. In 
turn, RIA neurons receive direct or indirect inputs from many 
sensory neurons, and are connected to numerous downstream 
interneurons, making them good candidates for transmitting 
information about ambient O2 to CO2 circuits. Previous work 
has shown that ttx-7 is required in the RIA neurons to promote 
appropriate synapse formation and to enable C. ekgans to navigate 
temperature gradients [38]. Expressing ttx-7 cDNA from the glr-3 
or gk-6 promoters, which drive expression exclusively in RIA [40], 
restored strong CO2 avoidance at 11% O2 (Figure 8D). By 
contrast, ttx-7 expression in AFD, using the gcjS promoter, or in 
AWC and AWB olfactory neurons, using the odr-1 promoter, did 
not. These data suggest that RIA interneurons are involved in 
communicating information from 02-sensing neurons and/or 
C02-responsive circuits, to enable its integration. 

We examined if CO2 elicited a Ca^""" response in RIA 
interneurons, and if this response was modulated by O2 context. 
We exposed animals expressing cameleon YC3.60 in RIA to a 
stimulus train in which we sequentially altered O2 and CO2 levels, 
and measured fluorescence changes in the cell body. 3% CO2 
evoked a Ca^"*" response in RIA neurons that was not significandy 
altered by background O2 (Figure 8E). These data suggest that 
RIA interneurons form part of a CO2 responsive circuit. Our 
inability to detect modulation of C02-evoked Ca^"^ responses in 
RIA by O2 levels could reflect a limitation of our imaging 



0 

-0.2 



B 



m -0.4-1 

I -o.eJ 



-0.8 
-1.0 



ID 

b 



0 i; 



-D -0.2 
c 

w -0.4 
X 

E 

(D -0.8 
5-1.0 



X 

a> 
c 

U5 
X 

S 
o 
E 

CD 



0 -t 



-0.2 
■0.4 
■0.6 



-0.8 



J 

NS 



■ 21 % Oj /I - 0% CO; ■ 21 % O; /3 - 0% COj ■ 21 

I I 11% 0,/1 -0% CO, □ 11% O,/3-0% CO, □ 11 



NS 

% 0, 15 - 
% O, /5 - 



0% COj 
0% CO, 



ttx-7: [promoter: :ttx-7cDNA] 




** 

■ 21%O,/1-0%CO, 
I I 11%O,/1-0%CO, 



5)^10 



t30 
•28: ° 

CO > 

-10 



|21%02 


21% O2 


21% O2 


11% O2 


11% O2 


11% O2 


21% 02 


0% CO; 


3% CO2 


0% CO2 


0% co_, 


3% CO2 


0% CO2 


0% CO; 






^f0m 








RIA 



0 



180 



360 



540 



720 



900 



1080 1260 



Figure 8. TTX-7 acts in RIA interneurons to promote CO2 
avoidance when ambient O2 levels are low. A-C. Mutations in ttx- 
7 strongly reduce CO2 avoidance at 11% O2 but have relatively weak 
effects on CO2 avoidance at 21% O2. ns, not significant, ** p<0.01. 
Student's f test. D. Expressing ffx-7 specifically in RIA neurons, using the 
glr-3 or glr-6 promoters, restores strong CO2 avoidance to ffx-7 mutants 
assayed at 1 1 % O2. Expressing ffx-7 specifically in AFD, using the gcy-8 
promoter, or in AWB and AWC, using the odr-l promoter does not 
rescue the ffx-7 phenotype. ns, not significant, ** p<0.01. Student's f 
test. E. CO2 evokes a Ca^* response in RIA neurons. Ca^* responses were 
measured in immobilized animals cultivated at 22°C using a pglr- 
6::YC3.60 Ca^* reporter. Shading highlights gas switch times. The CO2/ 
O2 stimulus train used is indicated above the plot. 
doi:1 0.1 371/journal.pgen.l 00401 1 .gOOB 

conditions. Alternatively, O2 could regulate RIA independently of 
Ca^""" entry, or could act on neurons downstream of RIA. 

Acclimation temperature and ambient O2 act 
combinatorially to regulate CO2 responsiveness 

Both acclimation temperature and acute ambient O2 concen- 
trations altered C. ekgans' responsiveness to CO2. We investigated 
how animals integrated information from all three homeostatic 
systems - temperature, O2 and CO2. We grew animals at either 
15°C or 22°C, and then assayed COj responses at 15°C or 22°C 
in the presence of either 21% or 11% O2. Our results suggest that 
the temperature and O2 sensing systems act additively to set CO2 
responsiveness. Decreasing O2 from 2 1 % to 11 % enhanced 
avoidance of 1 % CO2 regardless of acclimation temperature and 
assay temperature (Figure 9 A-C). Similarly, acclimating animals 
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to 15°C decreased avoidance of 1% CO2 at both 21% and 11% 
O2 (Figure 9A-C). As described previously (Figure lA), animals 
acclimated to 22°C avoided a l%-0% CO2 gradient more 
strongly when assayed at 15°C rather than 22°C. Changing O2 
from 21% to 11% further stimulated CO2 avoidance in these 
animals. These data highlight how C. elegans homeostatic responses 
are intertwined with each other. 

Discussion 

Previous acclimation temperature and current ambient O2 
levels set the aversiveness of CO2 to C. elegans. The temperature 
animals have experienced previously appears to modify CO2 
responsiveness by changing the CO2 receptive properties of AFD. 
Acute ambient O2 controls CO2 preference by regulating tonic 
signaling from the O2 sensing neuron URX. Changes in CO2 
responsiveness can be observed in shallow gradients with peak 
CO2 levels of 1%. Such gradients are likely to be ecologically 
relevant for C. elegans in the rotting fruit habitats where they are 
commonly found [41]. 
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Figure 9. Acclimation temperature and ambient O2 levels have 
additive effects on CO2 avoidance. A. Animals cultivated at 22"C 
but assayed at 15 'C avoid CO2 more strongly when ambient O2 is low. 
B-C. Reducing O2 levels from 21% to 11% increases CO2 avoidance 
regardless of acclimation temperature or assay temperature. In A-C, 
** p<0.01 , * p<0.05, ns, not significant, Student's f test. D. Coalitions of 
CO2 sensors elicit CO2 escape responses according to O2 environment, 
temperature experience, and CO2 stimulus dynamics. Triangles 
represent sensory neurons and hexagons interneurons. Black arrows 
indicate synapses. Several neurons respond to CO2 (blue arrows), each 
with distinct kinetics. Each of these neurons also responds to other 
sensory cues, as indicated. Three of the four identified CO2 sensors 
synapse directly onto the RIA interneuron. The fourth, AFD, synapses 
onto AIY which in turn synapses on RIA. The URX O2 sensor also 
synapses onto RIA. Note each neuron makes additional connections 
besides the ones highlighted here. 
doi:1 0.1 371/journal.pgen.1 00401 1 .g009 



C. elegans can thrive at temperatures that span ~1,5°C-25°C. 
Within this range, well-fed animals migrate to temperatures at 
which they were previously growing [13,29]. Temperature 
preference appears to be encoded in the AFD neurons: 
acclimation temperature changes the threshold at which rising 
temperature evokes Ca^"^ responses in this neuron [17,18]. We 
find that AFD neurons are required for temperature experience 
to change C. elegans' CO2 responsiveness. Acclimation temper- 
ature qualitatively reconfigures C02-evoked Ca^""" responses of 
AFD neurons. This re-configuration is retained in mutants 
defective in synaptic release, suggesting it can occur cell- 
autonomously. A speculative explanation of our observations is 
that AFD harbors multiple CO2 sensors whose contribution to 
the C02-evoked Ca^"*" response varies according to acclimation 
temperature. 

AFD neurons are exquisitely sensitive to CO2. They respond 
robustly to changes in CO2 that range from <0.01% C02/sec to 
>1% C02/sec. Remarkably, in animals acclimated to 22°C, the 
Ca^"^ responses evoked in AFD by slow (0.01% C02/second) and 
faster (0.04% C02/second) changes in CO2 are qualitatively 
dififerent. This may explain previous observations that AFD 
promotes CO2 avoidance in shallow CO2 gradients, but can 
inhibit CO2 avoidance in steep ones [28]. 

C. elegans avoid CO2 less strongly at high O2 than at low O2. 
Ambient O2 levels provide a contextual cue that modulates the 
aversiveness of CO2. We use the term 'contextual' because 
modulation can occur when O2 levels are constant, and is 
sustained over many minutes. Contextual modulation by O2 levels 
can be graded: as O2 decreases from 2 1 % to 1 1 %, CO2 avoidance 
rises. Modulation of CO2 avoidance by O2 requires the gcy-35 and 
gcj-36 soluble guanylate cyclases, which act in the O2 sensing 
neurons AQR, PQR and URX to transduce O2 levels, gcy-35 or 
gcy-36 mutants behave like animals kept at low O2, regardless of 
actual O2 levels. The activity of the URX neurons alone appears 
sufficient to inhibit CO2 avoidance at 2 1 % O2. Previous work has 
shown that URX neurons are tonically activated by high O2 [24], 
explaining the ability of these neurons to convey O2 context 
persistently to CO2 sensing circuits. 

Modulation of CO2 avoidance by O2 levels can be observed 
when N2 (Bristol), npr-1, glb-5(Haw); npr-1, or CB4856 (Haw) 
animals navigate l%-0% CO2 gradients. However, the degree of 
inhibition varies across these genotypes. In N2 animals, the 
inhibitory effect of O2 is limited by the action of the NPR- 1 2 1 5V 
isoform in 02-sensing neurons, npr-l 27 JFdoes not appear to alter 
the excitability of O2 sensors, since N2 and npr-1 mutants show 
similar 02-evoked Ca^"^ responses in URX, AQR or PQR ([22] 
and data not shown). Instead, we speculate that NPR-1 215V 
inhibits neurotransmission from URX, for example through G„ 
signahng [42,43], thus limiting the ability of URX to inhibit CO2 
responsiveness. Previous work has highlighted coupling of NPR-1 
215V to G„ in heterologous systems [44]. The potent O2- 
dependent inhibition of CO2 avoidance found in npr-1 mutants is 
suppressed by the gib- 5 (Haw) allele. This suppression appears to 
reflect a reduction in the excitability of URX. Tonic Ca^"*" levels in 
URX in glb-5; npr-1 animals kept at 2 1 % O2 was only as high as 
that found in npr-1 animals at 17% O2. In the CB4856 (Haw) 
strain the combination of the npr-1 215F and glb-5(Haw) alleles 
(potentially modified by other loci) enables a switch from 2 1 % to 
19% O2 to convert CO2 from a neutral to a strongly aversive 
stimulus. While this paper was in preparation independent work 
also highlighted modulation of CO2 avoidance by O2 in npr-1 
animals [45]. The assays used are different. Notably, in most of 
our work we used 1-0% CO2 gradients, whereas CarriUo et al. 
used 10%-0% gradients. 
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CO2 sensing in C. elegans is distributed across multiple sensory 
neurons, including the AFD and BAG neurons [28] (Figure 9D). 
Disrupting AFD and BAG abolishes CO2 avoidance at 21% O2, 
but CO2 avoidance at 1 1 % O2 is only partly reduced. Thus, CO2 
sensing neurons other than BAG and AFD can promote CO2 
avoidance at low O2. O2 modulation of CO2 responsiveness 
involves the RIA interneurons. ttx-7 mutants disrupt O2 modula- 
tion of CO2 responsiveness, and expressing ttx-7 cDNA selectively 
in RIA neurons rescues this phenotype. ttx-7 encodes myo-inositol 
monophosphatase. In ttx-7 mutants RIA neurons exhibit defects in 
localization of both pre- and post-synaptic components, including 
synaptobrevin, SYD-2 Liprin, and the glutamate receptor GLR- 1 
[38]. Synaptic communication via RIA is thus likely to be 
compromised in ttx-7 mutants, and may explain the O2/CO2 
integration phenotype. 

Previous studies of context-dependent changes in behavior in 
C. elegans have focused mainly on the effects of food or of food 
deprivation. C. elegans' migration in salt and odor gradients can 
switch from attraction to repulsion if animals are deprived of 
food in the presence of the chemical cue [46^9]. Food and food 
deprivation ha\'e also been shown to modulate C. elegans 
response to tcmpiuaturc' gradients [50]. It remains to be seen 
if acclimation temperature and ambient O2 levels have effects 
on other sensory modalities besides CO2 sensing. Whether CO2 
itself can act as a contextual cue regulating other C. elegans 
sensory responses, including thermotaxis and O2 sensing, is also 
unknown. 

The shallow CO2 gradients we study are likely to be common in 
the rotting fruit environments where C. elegans is frequendy 
found. However, the ubiquitous production of CO2 by 
aerobically respiring organisms means its value as a sensory 
cue likely depends crucially on context. Bacterial food, bacterial 
pathogens, predators, mates and conspecifics may all generate 
CO2 gradients. Context-dependence of CO2 responses has been 
observed previously. C. elegans CO2 responses are modulated by 
food, exposure to hypoxia, and starvation [25]. Moreover, not 
only context, but also the rate of change in CO2 concentration 
(whether it is slow or rapid), appears to modify the contribution 
of different C02-sensing neurons to C. elegans CO2 avoidance 
behaviors [28]. This complexity is mirrored in insects. For 
example in Drosophila airborne CO2 is aversive [51], whereas 
dissolved CO2 is attractive [52]. These properties are encoded 
by separate chemosensory neurons in the antenna (avoidance of 
gaseous CO2) and taste peg neurons (attraction to carbonation). 
Avoidance of airborne CO2 is inhibited by olfactory odors, 
presumably to enable flies to approach fermenting fruit [53]. 
Together, these data suggest CO2 sensing is remarkably 
sophisticated in both worms and flies. CO2 has been implicated 
in ageing in Drosophila [54], whereas 02-sensing neurons 
modulate longevity in Caenorhabditis [55], consistent with 
neurons sensing these gases also modulating physiology. 

Materials and Methods 

Strains 

Strains were maintained at 22°C with plentiful food using 
standard methods [56]. Strains used in this work are listed in 
Supplementary methods. 

Behavioral assays and analysis 

Spatial carbon dioxide gradient assays were performed as 
described, with slight modifications [25,28]. Briefly, rectangular 
PDMS chambers with a 33x15x0.2 mm space connected to gas 
syringes were placed over 100-200 worms on a 9 cm NGM agar 



plate. Assays ran for 20 minutes and the distribution of worms 
recorded by counting the number of animals in each of nine 
equal area divisions as well as in the two spaces at either end of 

the chamber. Animals were washed three times in a watch glass 
then transferred to the agar. A chemotaxis index was calculated 
by subtracting the number of animals in the low carbon dioxide 
half of the chamber from the number in the high carbon dioxide 
half and dividing by the total number of animals e.g. (A— B)/ 
(A+B), as shown in Figure lA. In chemotaxis assays, each data 
point represents the average of at least eight independent assays 
performed over three experimental days. Certified gases with 
indicated concentrations of O2 and CO2 were obtained from 
BOC UK Ltd. Assays marked 22°C were carried out at room 
temperature in a room in which temperature varied 22+/— 1°C. 
Assays marked 15°C were carried out in a small thermostat- 
controlled room set to 15°C. 

Statistical comparisons were carried out using the Student's t 
test or ANOVA, as indicated. 

Molecular biology and germline transformation 

Standard methods for molecular biology were used [57]. 
Cosmid and cDNA subcloning were performed using the 
Invitrogen Multisite Gateway Three-Fragment Vector Construc- 
tion Kit. 

Germline transformation was by microinjection [58] using 2— 
20 ng/|.il for the DNA to be tested, along with 50 ng/|j,l pJMZ-lin- 
15 (+) construct and carrier DNA, pBluescriptll SK (+). 

Ca^^ imaging 

Ca^"^ imaging was carried out as described previously [24,28], 
using an inverted microscope (Axiovert, Zeiss), a 40 x C- 
Apochromat lens, and MetaVue acquisition software (Molecular 
Devices). 

Supporting information 

Figure SI C02-evoked responses in AFD do not require the 
GCY-9 transmembrane guanylate cyclase (A, B), whereas BAG 
responses do (C, D). For all experiments animals were grown at 
22°C. 

(EPS) 

Figure S2 Disrupting gcy-33 reduces CO2 avoidance at 1 1 % O2, 
whereas disrupting gcy-32 or gcy-34 has no effect on CO2 
avoidance either at low or high O2. *, j()<0.05, **, ^<0.0I, ns, 
not significant. Student's Hest. 
(EPS) 

Figure S3 C02-evoked Ca^^ responses in ASE (A), BAG (B) and 
AFD (C) neurons are not altered by background O2 levels under 
our imaging conditions. CO2 and O2 stimuli are indicated above 
each plot. 
(EPS) 

Figure S4 ttx-7 mutants behave like N2 animals in 21%-0% O2 

gradients. 

(EPS) 

Text SI Strain list. 
(DOCX) 
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